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t^* ' Abstract 
ON 

The beta function of the vacuum energy density is analytically computed at the 
five-loop level in 0(iV)-symmetric 4 theory, using dimensional regularization in 
conjunction with the MS scheme. The result for the case of cubic anisotropy is also 
given. It is pointed out how to also obtain the beta function of the coupling and 
the gamma function of the mass from vacuum graphs. This method may be easier 
than traditional approaches. 



I Introduction 

In this article, we extend earlier work where the beta function f3 v of the vacuum energy den- 
sity in 0(iV)-symmetric 4 theory was computed at the four- loop level, to five loops. Integrals are 
dimensionally regulated and divergences removed by modified minimal subtraction (MS). For mo- 
tivations for this work, see We employ two different methods to arrive at our result for the 
0(iV)-symmetric case. In section |I], we use the scheme from [IJ] to determine the five-loop contri- 
butions to the vacuum energy renormalization constant Z v and to j3 v . All necessary subtractions 
result from consistency requirements while renormalizing the vacuum energy. As explained in 
this yields as a by-product 7 m through four loops and (3 g through three loops. In section |T|, we 
check the relevant recursion relations for the contributions to Z v . In section [VJ, we give f3 v through 
five loops for the 0(iV)-symmetric case. In section [V], we use a more traditional approach to check 
our result for (3 V . In section |VT], we extend it for the case of an additional cubic interaction. Our 
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work brings the evaluation of (3 V on a par with that of the other beta and gamma functions in 4 
theory (see @, |3| and references therein). 

For definitions and conventions, the reader is referred to the detailed article in [[!]. The only 
exception is (3 V itself, which we define here as 
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where the subscript B indicates that bare quantities are kept fixed. I.e. f3 v and (3 v ^ e in this article 
are just f3 v and /3„ j6 from Q divided by the renormalized coupling g. The connection to the constant 
Z v = 1 + h~ 4 J2T=i Zv,k{o) I ' e j which renormalizes the vacuum energy, is now 
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Our revised definition is more natural since even in the non-interacting theory, i.e. for g = 0, the 
vacuum energy density acquires a divergent part from the zero-point fluctuations of the uncoupled 
harmonic oscillators. It has to be renormalized and therefore runs with /i, such that 
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for the free theory. E.g., the renormalization group equation for the effective potential in d = 4 — e 
dimensions for h = reads now [ITT| 
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For the free theory, we have (3 g ^ = 7 mj£ = 7^ e 
energy 

V(g = O,m 2 ,0 2 = 0,h = 0,/i 2 
with the MS renormalization scale 
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and, going to four dimensions, the MS vacuum 
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f = 4vr / u 2 e- 7B (6) 
is immediately seen to satisfy (|J) with (3 V given by (|]) with no spurious appearances of a coupling. 



II Five-Loop Contribution to Z v and f3 v 

In this section, we determine the five-loop contributions to Z v and to j3 v in the 0(iV)-symmetric 
model. We proceed in the same way as in |1|]: We employ modified Feynman rules, where the 
one-loop mass correction is absorbed into a modified mass in the propagator. This reduces the 
full set of diagrams in table [I] to the reduced set in table 0, where we have also included the six- 
loop diagrams, whose evaluation we save for a later day. We arrange that the appearances of the 
wave function renormalization constant in the propagator and coupling cancel from the outset. 
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Table 1: Vacuum diagrams through five loops and their symmetry factors. In the equations in the 
text the symmetry factor is considered part of each respective diagram. 



There are no extra counterterm rules, since all counterterms are already contained in the Feynman 
rules for vacuum energy, propagator and coupling; i.e. we use bare values for vacuum energy, mass 
and coupling in the Feynman rules and expand our results for diagrams to the necessary order in 
the renormalized coupling g. For this purpose, we have to formally construct the renormalization 
constants Z m , and Z g , for mass, modified mass and coupling, through four, four and three loops, 
respectively, from 7 m and (3 g [0. 

The calculation through four loops is detailed in PJ . Continuing to five loops, we keep all divergent 
terms in zero through five loops through order g 4 . The relevant five-loop diagrams are 
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Table 2: Remaining diagrams through six-loop order after absorption of the one- loop mass correction 
into a modified mass in the propagator, i.e. after a careful resummation of the quadratic part of the 
Lagrange density. 
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where Is a , 1^ and are defined, their evaluation sketched and their results given in appendix [A[ 
Using these, the formally reconstructed Z m , Z m and Z g , as well as the results of one gets 
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Demanding the cancellation of logarithmic terms gives the three-loop coefficient of f3 g and the four- 
loop coefficient of j m , 

a 33JV 2 +922iV+2960+96(5A r +22)C(3) /,-,\ 

P3 - 432 , (J-J-J 

(Af+2)W 2 -7578Ar-31060-48(3Ar 2 +10Ar+68)C(3)-288(5iV+22)C(4)] /-, n \ 

«4 - 15552 , (42) 

which coincide with known results (see, e.g., 0, ||). Demanding subsequently ( |l0l) to be finite as 
e — > gives 
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The relation (fj) then gives us the five- loop coefficient of (3 V = J2kLi ^fc5 ,fc_1 through 5 5 = 5Z" 5 /2. 



Ill Recursion Relations for the 

The relevant recursion relations among the renormalization constants are M 
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Using the results from H, it is straightforward to check that all of the above relations hold. 
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IV (3 V Through Five Loops 



Combining our result for 65 with the four- loop result from Jj]], we arrive at the five-loop result for 
the 0(iV)-symmetric case, 

R ( n \ - N I N ( N + 2 ) n 2 I ^(^+2)(iV+8)[12C(3)-25] 3 
HvKy) — 4 "r 96 » "r" 1296 » 

, jV(jV+2)[-319jV 2 +13968Af+64864+16(3Af 2 -382jV-1700)C(3)+96(4Af 2 +39jV+146)C(4)-1024(5Af+22)C(5)] 4 
~T~ 82944 » 

+0(^ 5 ) . (22) 

We note that among the beta and gamma functions, j3 v is the easiest to compute. Therefore, 
it would be the prime candidate for the first complete non-trivial six-loop calculation, since only 
the last six diagrams in table |] have to be computed, of which the first three are essentially trivial. 
After converting the quartically divergent integrals into logarithmically divergent ones by twice dif- 
ferentiating with respect to m 2 (see appendix |A|), a subset of the diagrams necessary for the six-loop 
renormalization of the coupling has to be evaluated. Since this subset consists of the diagrams where 
the four external lines are attached to only two vertices, the most difficult topologies are absent and 
the computation should be considerably more easy than the full coupling renormalization at this 
level. A similar statement is true for the comparison with mass and wave function renormalization. 

We further note that the fc-loop coefficients of 7 m we get as by-products here and in [Q 
come from e _1 ln(m 2 //i 2 ) terms — and therefore originally from e -2 terms — of (k + l)-loop vacuum 
diagrams. Similarly, the fc-loop coefficients (3k of (3 g we get come from e _1 ln 2 (m 2 //i 2 ) terms — and 
therefore originally from e~ 3 terms — of (k + 2)-loop vacuum diagrams. These may be easier to 
compute than the e _1 terms of generic /c-loop diagrams contributing to (3% and a^. In particular, and 
in contrast to 0, we do not have to consider the topology of eight-loop vacuum diagrams to obtain 
the five-loop contribution to (3 g , but only to compute the 1/e 3 terms of seven-loop vacuum diagrams. 
To obtain the five-loop contribution to j m , we have to to compute the 1/e 2 terms of six-loop vacuum 
diagrams. 

However, in the current approach of evaluating loop integrals (to be described in appendix 0), the 
divergences related to the two-point function at one lower loop (and therefore contributing to 7 m at 
one lower loop) and to the four-point function at two lower loops (and therefore contributing to (3 g at 
two lower loops) than the considered vacuum graph integrals need to be evaluated as sub divergences 
of these vacuum graphs anyway. Therefore, when evaluating the integrals as described in appendix [A], 
no advantage has been gained for computing 7 m and j3 g , at least as far as the necessary computational 
techniques are concerned. 



V (3 V and Z v from the Standard Method 

Along more traditional lines, we have 

Z v — 1 + ^4 KLRYj vacuum graphs , (23) 

where R removes sub divergences of graphs and K, isolates the negative powers in e (see for a 
clear definition of /C and R). The sum in ( P3"D goes over the graphs in table |TJ and the Feynman rules 
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to be used are 



a — b = -^- 2 , (24) 

pZ _|_ m z 

a \/b _ (4vr) 2 ^ 

X , = -[tiabdcd + Sachd + 8 ad 5 bc ] — - — (25) 
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with renormalized quantities m 2 and g, since all sub divergences are removed by the R operation. We 
have carried out this program from one through five loops in order to also check the results in [p]]. 
Using the results of appendix [FJ, it is easy to see that we get the same Z v and therefore also (3 V as in 
the approach above. 

If a diagram D separates into n diagrams D k with independent integrations, then the standard 
definition of ICR yields 

n n 

1CRD = ICR J] D k = (-l) n+1 [] KRD k ■ (26) 
fc=i fc=i 

Since /C picks out only the pole terms in e, the 1/e terms of Z v and therefore also (3 V receive only 
contributions from diagrams whose integrations are not separable. Thus, for the calculation of f3 v , 
it is again sufficient to consider the diagrams of table 0. The two-loop diagram is no longer needed, 
though. However, since, within this program, separable diagrams can be simply put algebraically 
together from lower-loop diagrams, including them in the calculation provides together with the 
recursion relations between the Z x 111] a convenient cross-check in the determination of Z v and f3 v . 



VI Cubic Anisotropy 

Here we give /3 V for the case of a cubic anisotropy |3], |§. We rename g — > g\ and introduce a second 
coupling g 2 through 

£ = ^d^Bid^Bi + l -m\$Bt§Bi + (gisT^li + 92 B T^ fa^kfa + > ( 2T ) 

where repeated indices are summed over (/i = 1, . . . , d and i,j,k,l — l,...,N), the subscript B refers 
to bare quantities and where 

T ijli = g(^j4z + S ik 8ji + diidjk) , (28) 

T (2) _ x _ f 1, i = j = k = l, . v 

T m - 5 ijkl = j Qj elge _ (29) 

Using standard methods [|J, one arrives at 

i / dz Vjl dz v A 

Pv = — \9i^ — + 92^ — • (30) 



2^i V 9 9i dg : 
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The seemingly asymmetric treatment of gi and gi is entirely due to our definition of the constant 
term in (|2"7j), since the definitions (^Sj) and (^) are not used for the derivation of fl3"0|). However, 
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it is easy to see that e.g. replacing the constant term in ( |2T| ) by m%IiB / '[(47r) 2 g2 B ] and redefining Z v 
accordingly exchanges gi and #2 in (T30|) , but leaves /3 V invariant. 
Through five loops, we get 
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2304 9l92 
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9216 ^2 1 
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which reduces to (UT) for , c/2) (fisO). 



Note Added 

Recently, the result fl22|) has been confirmed in an independent calculation by S.A. Larin, M. Monnig- 
mann, M. Strosser and V. Dohm, |cond-mat/9711069| , where also its application to three-dimensional 
systems is considered. 



Acknowledgements 

I am grateful to Stefan Bornholdt for numerous helpful discussions. This work was supported by the 
Deutsche Forschungsgemeinschaft (DFG). 



Appendix 

A Five-Loop Integrals 

Our strategy for computing the five-loop integrals I^ a , i^, he is different from the way of computing 
integrals in Here we take two derivatives of each integral with respect to m 2 to convert it 
into a sum of logarithmically divergent integrals. Then we subtract sub divergences using the ICR 
operation, so that we can set most masses in propagators in the resulting expression to zero without 
changing its divergent part. This allows us to do the relevant integrals using the methods of infrared 
rearrangement ||, (our own modified version of) the R* operation and the results of @, being 
in turn partially based also on the integration-by-parts algorithm ||. For an introduction to these 
techniques and the multi-loop renormalization of 4 theory in general, see 0. At last, we evaluate 
the terms containing sub divergences that we have subtracted above and add them again. These 
terms have less than five loops and can be computed by either recursively continuing this procedure 
or by using the results from 
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One might argue that this method of computing integrals is not independent of the standard 
method referred to in section |V|. Nevertheless, it provides us with many more cross-checks as de- 
scribed in sections [TT] and |TT[ 

All diagrams in this section only refer to momentum space integrals without symmetry and group 
factors. To get the notation in line with appendix [FJ, we rename the integrals /f and J| c from |TJ 
I 3a and J 4a , respectively, so that their finite parts defined in |J are now I 3a j = I%j and J 4a j = if j. 

A.l I 5a 

Define 7 5a by 
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kpqrs 



(p 2 +m 2 )[(fc+p) 2 +m 2 ] (q 2 +m 2 )[(k+q) 2 +m 2 ] (r 2 +m 2 )[(k+r) 2 +m 2 ] (s 2 +m 2 )[(k+s) 2 +m 2 ] 
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Since I 5a oc (m ) 2 



(m 2 ) 2 2<E we can write 
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The three integrals on the right hand side can now be evaluated as described above. The result is 
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120 In 2 s£ - 424 In ^ + + 24((2) 
116 In 3 s£ - 298 In 2 + (258 + 156C(2)) In s£ + (- ^ - 170C(2) + f C(3) 
f In 4 - 45 In 3 s£ + (i§2 + 45C(2)) In 2 fr - (f + 71C(2) + 12C(3)) In ^ 
+ (-|f + f C(2) + f C(3) + ^C(4) + f C(2) 2 )] } 

+ {4 ^4 e2 ha,f + J^yrJiaJ + ha J , (34) 

where I 5aJ = 0(e°). 



A.2 J 5b 

Define /5b by 
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(fc 2 +m 2 ) 2 [(fc+p+g) 2 +m 2 ](p 2 +m 2 )(g 2 +m 2 )[(A:+r+s) 2 +m 2 ](r 2 +m 2 )(s 2 +m 2 ) 



(35) 
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We can write 
hb 
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and evaluate this to give 
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^10 | 5e 5 "I" e 4 
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mf, 452 
/I 2 5 



+ i 



+ - 



+ 
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where J 5bi/ = O(e ) 

A.3 I 5c 

Define I§ c by 



120 In 2 - 226 In + + 24((2) 
46 In 3 f + ifi In 2 f + (if - 6C(2)) In f + (-^ - f C(2) + 4C(3))] 
-f ln ' W + f ln3 # " + f C(2)) In 2 f + + ?C(2) - 10C(3)) In f 
+ - f C(2) + f C(3) - |C(4) + |C(2) 2 ) 
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and evaluate this to give 
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^ = (^{0 + i[-481n^ + fl+^[6Oln 2 ^- 2281n^ + ^ + 12C(2) + fC(3)] 



58 In 3 s£ - 129 In 2 + + 78((2) - 16C(3)) In § 
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+ (^-81C(2) + fC(3) + fC(4) 
+1 [f In 4 f - f In 3 f + (JL + f C ( 2 ) + 20C(3)) In 2 f 
+ (-1 - f C(2) - 26C(3) - 12C(4)) In f 

+ f C(2) + |C(3) + fC(4) - f C(5) + f C(2) 2 + 4C(2)C(3)) ] } 



(40) 



where L, 



5c,/ 



B /C.R on Quartically Divergent Diagrams 

Using standard methods, we have evaluated the quartically divergent diagrams from table [I] with 
sub divergences removed, as needed for the determination of Z v in section [V]. As in the last section, 
the diagrams in the following refer to momentum space integrals, while symmetry factors S x and 
group factors G x (gi,g 2 ) are written out separately. As indicated, we define the group factors to 
contain not only the contributions from the contraction of the 8y and Sijki tensors, but also their 
accompanying factors —gi/3 and —g%, respectively. For diagrams with separable integrations, we only 
give the 0(iV)-symmetric result G x (g,0), since these diagrams are exclusively used for cross-checks, 
which we have performed only for the 0(iV)-symmetric case. For diagrams whose integrations do 
not separate, we give G x {g\,g2). 

Our strategy to compute the quartically divergent integrals with subdivergences removed is to 
convert quartically divergent integrals by twice differentiating with respect to m 2 to logarithmically 
divergent ones and subsequently evaluate these by the methods of 0-|§. We exploit the fact that the 
operator d/dm 2 commutes with ICR (in the same way as differentiation with respect to an external 
momentum commutes with ICR, see |]TU[) and that the result of ICR acting on a quartically divergent 
diagram is proportional to m 4 , 

We have evaluated all logarithmically divergent diagrams encountered en route and found agreement 
with the respective results in 0. Therefore, we do not display the results for the logarithmically 
divergent diagrams here. In order to also check our results in we have carried out this program 
from one through five loops and give the results for all diagrams of table [I] below. 

The notation of I x refers to the integral I x with subdivergences removed by ICR. The integrals 
I 3a and J 4a below are identical to Iff and If in M, respectively. 




B.l One Loop 

I la = ICR O = W^R 0=wk' ( 42 ) 
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5ia — \ i 
G la ( gi ,g 2 ) = N. 



(43) 
(44) 



B.2 Two Loops 



I 2a ee ICR GO = m'JCRQQ = 



5*2a — | 



G 2o (<7,0) = -±JV(iV + 2)<7. 



4,2 5 



(45) 

(46) 
(47) 



B.3 Three Loops 



ha = ^C-R 1 



a _ J_ 



G 3a ( gi ,g 2 ) = $N(N + 2)g\ + 2N 9l g 2 + Ng : 



(48) 

(49) 
(50) 



$3b = ^ , 

G 36 (<?, 0) = lN(N + 2) V ■ 



8m 4 
( 47r )6 e 3 



B.4 Four Loops 



(51) 

(52) 
(53) 



1 4a 



ICR 



(47r) fc 




\m 4 lCR 
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+ 24 




(54) 



a _ J_ 

*J4a — 4§ j 



G4o(^i, 52) = -^iV(iV + 2){N + - |iV(iV + 8)^ 2 - 3iV M 2 - Ng* 



(55) 
(56) 
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Sib — , 

G 46 (ff,0) = -iiV(JV + 2)V- 



(58) 
(59) 



16m 



(47r)8e 4 



C _ J_ 

° 4c ~~ 32 ' 

G 4c (g,0) = -±N(N + 2) 3 g 3 - 



(60) 

(61) 
(62) 



B.5 Five Loops 




0, 



5. 



■id 



48 



G 4d {g,0) = -^N(N + 2) 3 g 3 ,- 
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(64) 
(65) 



^5a = ^-R 




1^4 



rrfKR 
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(47T) 1 



192 352 , 144 , 1 f 168 272 /-(o\\ , 1 ( _ 319 , 12 a/q\ , 96a//1n\ 
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G 5a (9i,92) = ±N(N + 2)(N 2 + QN+20)gt + ±N(N 2 + QN+20)g!g 2 + l:N(N 
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+ 24 ■ 



+ 6 




m 4 ^192_208_J4__|_138_ 71 \ 

(47r) 10 V5e 5 5e 4 5e 3 5e 2 20e J ' 



(69) 



c _ l 

°5b — 144 



(70) 

G 5b ( gi ,92) = lN(N + 2fg\ + |iV(iV + 2)g\g 2 + |iV(iV + K)g\g\ + 4iV^ 2 3 + A^ 2 4 . (71) 



13 



he = K<R 



= \m A K,R 
















! +4. 






b 




















+8. 


be 






D 


3?. 















(47T) 1 



96 _ 304 , 1 ( 1592 , 192 a/o\\ , 1 / 340 , 24a/o\ V§/-(A\ 
5e 5 5e 4 + e 3 ^ 15 + 5 CAOJJ -t- e2 ^ 3 + 5 5 

+i - 16C(3) + f C(4) - f C(5)) 

c _ J_ 

° 5c ~ 32 ' 



(72) 



(73) 



G 5c (9i,92) = ±N(N + 2)(5N + 22)gi + ±N(5N + 22)g!g 2 + lN(N + 26)g 2 1 g 2 2 +Wg 1 gl + Ngl. (74) 



l5d 



ICR 



(47T) 1 




\m 4 lCR 




+ 2 




+ 



S-g + I + M~f + 4 C( 3 ))]> 



(75) 



G 5d (g : 0) = ±N(N + 2y(N + 8)g±. 



(76) 
(77) 



( 5e 



ICR 




m A KR 




(47r) 10 V3e 4 e» 



C _ J_ 

° 5e ~~ 48 ' 

G 5e (g,0) = ±N(N + 2)Y- 



(78) 

(79) 
(80) 



h f = ICR 




= m 4 )CR 




m 4 / _32_ 32 , 8 

(4^p V3e 5 3? + P 



(81) 
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s 5f - h ' 

G 5f (g,0) = ±N(N + 2)Y 



(82) 
(83) 



\m 4 lCR 





(4 



to 4 / 32 80 I 2 
tt) 10 3e 4 ^ e 3 



<? - -L 

5 9 48 ' 



G 6ff O/,0) = iJV(JV + 2)V- 



(84) 



(85) 
(86) 



CXXXX) =- 4 ^CXXXX) 

GBfc(»,0) = iJV(N + 2)V- 



32m 4 
( 47r )iO e 5 , 



(87) 

(88) 
(89) 



hi = ICR 




0, 



C _ 1 
°5i — 32 , 



G u (g,0) = ±N(N + 2)Y 



(90) 

(91) 
(92) 




0, 



C _ 1 

128 ' 



G Sj (g,0) = ±N(N + 2)Y 



(93) 

(94) 
(95) 
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